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ABSTRACT 

We present C-functions for static and spherically symmetric spacetimes in Lovelock gravity 
theories. These functions are monotonically increasing functions of the outward radial coor- 
dinate and acquire their minima when evaluated on the horizon. Unlike the case of Einstein 
gravity, where there is a single C-function, we find that this function is non-unique in the 
case of Lovelock gravity. We define two C-functions, which agree at the horizon giving the 
black hole entropy, and state the different energy conditions that must hold in order for these 
functions to satisfy the monotonicity condition. 
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1 Introduction 



Since the remarkable discovery [1] that a black hole has entropy proportional to the area of 
the horizon 

many approaches have been proposed to count the number of quantum mechanical states 
that contribute to this entropy. Particularly intriguing are connections with two dimen- 
sional conformal field theory. In this context Solodukhin [2] and Carlip [3] showed that 
if we consider the black hole horizon as a boundary condition on the radial fluctuations 
of the metric then we obtain, in the vicinity of the horizon, an infinite-dimensional group 
of conformal transformations in two dimensions with corresponding Virasoro algebra that 
contains the Bekenstein-Hawking entropy as a central charge. For quantum theories in two 
dimensions, Zamolodchikov [4] was able to prove a set of properties satisfied by what is 
called the C-function under renormalization group flow. This function was shown to be a 
function of the couplings of the theory that is monotonically decreasing as one flows to lower 
energies. For fixed points of the fiow, corresponding to the extrema of this function, the 
C-function reduces to the central charge of the Virasoro algebra. In [5] a holographic version 
of Zamolodchikov's C-theorem was proved by studying the renormalization group fiow along 
null geodesic congruences in asymptotically AdS spaces. Further, Sahakian [6] proposed a 
covariant geometrical expression for the C-function for theories which admit a dual gravita- 
tional description. In this description, the IR region is deep interior and the fiow outward 
in radius is toward the UV region in the QFT sense. Another possible interpretation of the 
holographic picture using the moduli fiow, in the context of the attractor mechanism, was 
given in [7]. 

It was shown by Goldstein et al [8] that in 4 dimensional Einstein gravity, coupled to 
matter fields that satisfy the null energy condition, one can define a simple C-function for 
static asymptotically fiat solutions. This function is given by 

c>,^^, (.) 

where A{r) is the area of the two sphere as a function of the radial coordinate and G4 is 
Newton's constant in 4-D. It was proved in [8] that the equations of motion imply that A{r) 
must decrease as one moves inwards from asymptotic infinity. Also, C(r) coincides with the 

entropy at the horizon. 

On the other hand the entropy in higher curvature gravity is given by the integral of a 
particular local quantity on a spatial cross section E of the event horizon [9, 10] 

S ^-2n / — eabecd\^dn , (3) 

OKabcd 

where L is the Lagrangian, Sab denotes the binormal to the horizon cross section and Vh dQ 
is the volume element induced on E. 

The question of whether one can define analogous C-functions in higher curvature gravity 
was raised in [11], where it was shown that a similar C-function can be obtained by evaluating 
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Wald's expression for the entropy (3) on a general spacelike surface instead of a spatial cross- 
section of the event horizon. Although the authors in [11] pointed out the monotonicity of 
the C-function for f{R) gravity (this had been established earlier in [12]), they were not able 
to draw a conclusion about the monotonicity of such functions in a general theory of gravity. 

An important class of higher curvature gravity theories is known as Lovelock gravity [13]. 
These are the most general second order gravity theories in higher dimensional spacetimes. 
A general formula for the entropy of stationary black holes in these theories was obtained by 
Jacobson and Myers [14]. It includes a sum of intrinsic curvature invariants integrated over 
a cross section of the horizon. This entropy coincides with the result one obtains using the 
Wald's formula ^ (3). In addition, similar calculations to [2] were performed for Lovelock 
gravity [17] and it was shown, as in [2] for Einstein gravity, that the central charge of the 
Virasoro algebra is proportional to Jacobson-Myers entropy. 

One can ask whether a C-function similar to that of [8] exists for the static spherically 
symmetric black holes in Lovelock gravity. In other words, one asks if a monotonically 
increasing function of the outward radial coordinate may exist under certain conditions, 
which reduces to the entropy when evaluated on the event horizon. We address this question 
in the present work. We show not only that such a C-function exists, but also that this 
function is non-unique. In fact we find two different C-functions that we call C-functions 
of the first and second kind. These functions exist provided that the matter content satisfies 
respectively the null, as in [8], and the weak energy condition, and that the spacetime is 
asymptotically flat. 

In the next section we review the construction of pure Lovelock gravity. Then we re- 
view an argument proving the monotonicity of A{r), the area of concentric spheres, in the 
spherically asymptotically fiat spacetimes. In section 3 we introduce the C-functions of the 
first and second kind of pure Lovelock gravity theories and we prove the monotonic behav- 
ior of these functions. Then in section 4 we consider the case of general Lovelock gravity 
and the behavior of the general C-functions in this theory. The proof of the monotonicity 
for C-functions of the second kind in general Lovelock gravity is cumbersome and requires 
thorough analysis for general polynomials of arbitrary degree. Wc work out the proof for 
Gauss- Bonnet gravity and we present numerical results for the second and third order Love- 
lock theories. These numerical results confirm our analytical results in the Gauss-Bonnet 
case and indicate that the C-function of the second kind may be monotonic in a general 
third order Lovelock gravity theory as well. 



2 Lovelock gravity 

The Lagrangian density for general Lovelock gravity in D dimesnions is £ = Z^j^io (^m^m, 
where Cm is given by 

r ^ /~r raifel-.-ambm p Cldl p Cmdm (A\ 

~ 2"* y "ci<ii...c^d^ -"-aioi ■■■■■^^-ambm ' 

am is the m'th order couphng constant, [-D/2] denotes the integer value of D/2 and the latin 
indices a,b,c and d go from to D — 1. The 5 symbol is a totally antisymmetric product 

■^The black hole entropy for Lovelock AdS gravity can also be obtained directly from a background- 
independent regularization of the Euclidean action as was shown in [15, 16]. 
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of 2m Kronecker deltas normalized to take the values of ±1. The term Cq — is the 

cosmological term, while L\ — \J—g 5"^^^^ i?ai6i'^^''V2 is the Einstein term. In general is 
the Euler class of a 2m dimensional manifold. 

As a special class of general Lovelock gravity, we take a theory with highest order in- 
teraction m < [{D — l)/2], and send the coefficients of all the lower order terms to 
zero. We call these pure Lovelock gravity theories, with pure Einstein gravity as the first 
non-trivial example [18]. There has been intensive effort to study black holes as well as their 
thermodynamic properties in the context of Lovelock gravity ( see e.g. [19]- [27]). 

In the following we will be interested in static spherically symmetric spacetimes. Hence, 
the metric can be assumed to take the form 

ds^ = -a(r)2 dt'' + ^ + b{rf dQl , (5) 

where rffi^^ is the metric on the unit n = D — 2 sphere. The nonzero components of the 
Riemann tensor for the above metric are given by 

it:,/* = -ia"a + a"), R,f^ = -a'{^-^ + ^)6i, (6) 

ab 

Rti'^ = -"^Si, Rr,''-^-^p-5^. (7) 

The equations of motion following from the Lovelock Lagrangian (4) are given by Qj ^"^^ 

_2m+l rjnf ^J^g^g 

r; f (m) _ i-faibl...ambm td cidi p Cmdm / o\ 

and T/ is the energy-momentum tensor. Using eqs. (6) and (7), and through repeated 
application of the identity 



^t-X = 2{D-p + l){D-p + 2) CX-t , (9) 



we obtain 



^,(^) _ 2"^+^mn! /b" ^ a'b'\ f l-aH'^ y-' ^ 2-n! fl-a^f/^Y 



-yr (m) 



(n-2m + l)! \b ab J \ 6^ j {n-2m)\\ 6^ j ' 

^7Z-^A^J^\ ■ (10) 



(n-2m + l)! V / V / (n-2m)!V 



2.1 Monotonicity of h and C-function in Einstein gravity 

As a warm up, we recall how the results of Goldstein et al [8] come about. For m = 1 we 
obtain the special case of Einstein gravity in D dimensions for which eqs. (10) take the 
simple form 

^/(^) = -4-^(^ + ^J+2n(n-l)^^^j, 

^/(^) = -4nM+2n(n-l)(l^). (11) 
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Now, consider the coupling of this theory to matter that satisfies the null energy condition, 
i.e. the energy-momentum tensor satisfies the condition 

Tabee>o (12) 

for all null vectors As a special case, one may take a perfect fluid with energy-momentum 
tensor given by 

Tab^{p + p)UaUi,+pgab, (13) 

where p, p and [/„ are respectively the fluid energy density, pressure and D-velocity. In this 
case the condition (12) reads p + p> 0. Goldstein et al [8] showed that in this system 6 is a 
monotonically increasing function of r for any static, spherically symmetric and asymptot- 
ically flat spacetime. To show this one takes a particular linear combination of Qt^^^^ and 

Gt' - G/ = -64 TT Gd {Tt' - T;) = 64 tt T^b CC' > , (14) 

where ^" = {CO are components of a null vector, satisfying the relation, (^*)^ = — gf** and 
iO"^ — 9^^ 1 ^^'i have used ai — I/IQttGd and Go is Newton's constant in D dimensions. 
Using eqs. (11) one obtains 

a'b"^-'-^^bT^,ee. (15) 
n 

As long as we are outside the horizon, > 0, we obtain b" < 0. Given that the spacetime 
is asymptotically flat, Goldstein et al [8] were then able to show that b{r) is a monotonically 
increasing function of r. The additional steps required are given below in the context of pure 
Lovelock gravity. 

Using this fact, we see that a possible C-function for Einstein gravity takes the simple 
form 

= h-^l^ Gd = Asn/^ Gd , (16) 

where is the volume of the unit n sphere. This expression coincides with the entropy 
when evaluated on the horizon. 



2.2 Monotonicity of h in pure Lovelock gravity 
The null energy condition 

Consider pure Lovelock gravity in D dimensions and of order m < [{D — l)/2] coupled to 
matter that satisfies the null energy condition. We show that in such a system 6 is a mono- 
tonically increasing function of r for any static, spherically symmetric and asymptotically 
flat spacetime given the positivity of the coupling constant, i.e. > 0. As we did in the 
case of Einstein's gravity, we take linear combinations of Gt' ^'^^ and G/ ^"^^ 

Gt' ^"^^ - G/ ^"^^ = (r** - r;) = — Tab re' > . (17) 



Using eqs. (10) we obtain 



^■'^-^^.^^J^r (^T.u'e). (18) 
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Now, let us assume that the metric (5) describes a black hole. We assume cosmic cen- 
sorship, so that 6(r) 7^ on, or outside the horizon. Without loss of generality, we can 
then assume that 6 > on the horizon. Asymptotic flatness is consistent with b{r) ^ ±r 
as r — > 00. However, the case b{r) ^ — r depicted on the right in figure (1) is ruled out by 
our assumption of cosmic censorship. The other two graphs in figure (1) make clear that for 
b{r) to then fail to be monotonic between the horizon and infinity, it must have at least one 
minimum in this range. However, this is ruled out by eq. (18). Assume 6' = at some radius 
ro, since > outside the horizon, it follows from eq. (18) that b"{ro) < 0. Therefore, local 
minima are not allowed. This proves the monotonicity of b for pure Lovelock gravity. 




Figure 1: Different possibilities for the function b{r). 



The weak energy condition 

The weak energy condition states that the energy density of any matter distribution, as 
measured by any observer in spacetime, must be nonnegative, i.e. 

TabX''x'>0 (19) 

for any future-directed timelike vector x". For our static, spherically symmetric spacetimes, 
we can take x° fo have only one non-vanishing component x* = 1- Hence, using (19) we 
obtain = —a?T^ = p > 0. Moreover, using the energy-momentum tensor of a perfect 
fluid (13) we find p + p > 0. This shows that the null energy condition (12) is a special case 
of the weak energy condition. 

Adopting the weak energy condition, we can show below that b" < 0, not only at the 
local extrema, but for all r where th < t < 00. To this end we rewrite the first equation in 
(10) in the form 

^{b''Tm,n) = -bWI, (20) 

dr 

where 

r = ^ Ji^^r (21) 

2(n-2m + l)\ [ b^ i ' 



Integrating eq. (20) and using T/ = —p we obtain 

' \ 1+71— 2m „/ „ 

b{r) J a-m n\ 6(r)^+"' 



(1 



/ \ 1+71— 2m „/ „ N , 
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We have shown above that the null energy condition ensures that 6(r) is a monotonic function 
and hence b'{r) > for th < r < oo, given that b(r) > 0. Moreover, using the weak energy 
condition, p > 0, we see immediately that the l.h.s of eq. (22) above is always greater than 
zero. It is trivial to see that the same result holds for b{r) < 0. Finally using (18) proves 
that b"{r) < for all r, th < r < oo. This is analogous to the result of pure Einstein gravity 
eq. (15) although in the latter case one uses only the null energy condition. 



3 Entropy of pure Lovelock black holes and C-functions 

The C-functions we want to find should by definition reduce to the black hole entropy when 
evaluated at the horizon. We begin this section by recalling the expression for the entropy 
in Lovelock gravity. 

A general formula for the entropy of stationary black holes in Lovelock gravity was 
obtained by Jacobson and Myers [14] using Hamiltonian methods. They showed that the 
entropy of a black hole in pure Lovelock gravity of order m is given by 

= An mam I d"x V^Cm-i{h) , (23) 

where h is the induced metric on the horizon, and the integration is evaluated on any 
spacelike slice of the Killing horizon. Using the spherical components of the Riemann tensor 
Rij^\h)^ calculated from the induced metric on the horizon, and the basic definition of the 
Lagrangian in eq. (4) along with the second eq. in (7) and the identity (9), we obtain 

Qim) ^ 4 7rmn!a^ , n-2m+2 

(n-2m + 2)! ^ ■ K ) 

Note that for m = 1 with ai = l/167rG4 we obtain = f2„6"/4G' = A/4G'4, where A is 
the surface area of the horizon, which is the Bekenstein-Hawking expression for the entropy 
of Einstein gravity. 

One can also use Wald's expression for the entropy (3) to obtain the same result (24) 
above. To show this we start from the Lovelock Lagrangian L^, where Cm — ^/—gLm and 
Cm is given by (4), to obtain by direct calculations 



dL _ 2m 

gj^ /1/2 ~ 2m '^Cldl...Cm-ldm-lhh 



rai6i...a„i_i6m_ieie2 p cidi / N p Cm-idm-1 I „\ ( 

"cAdu..Crr,.-\drr,.-\hfi ^a,ibi ) ■ ■ ■ ■ ^am-ibm-1 yy)- l^^J 



Using Ert — 1, the only nonvanishing component of the binormal to a spacial two-surface 
concentric with the horizon, we find 

r)T 9 ni 

f: Jlh — rnjl-^m-ljm-l p. . klh ( \ p. . fc^_lZ^_l/ \ 

^"^f-^ — Dm "klh...km-llm-l ^■^in •••• -"■lm-Um-1 {9)^ l^^J 



where Rij (g) denotes the spherical components of Riemann tensor calculated from the full 
metric. Using the last expression together with eq. (7) in Wald's formula we obtain 

I / -1 2 y2 \ "1—1 
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where b and the bracket are to be evaluated on the horizon, i.e. at a = 0, and hence we 
reproduce the same result given by eq. (24). 

A C-function should extend the entropy of spheres on a constant time slice away from 
the horizon. We see that eqs. (24) and (27) suggest different ways of doing this. This lead to 
two possibihties for C-f unctions. What we call the C-function of the first kind below is based 
on eq. (24), while the C-function of the second kind is based on eq. (27). The difference is 
that eq. (24) involves only the intrinsic curvature of the spheres, while eq. (27) includes the 
extrinsic curvature as well. 

3.1 C- functions of the first kind 

In [8] it was shown that one can take the C-function in Einstein gravity, for the static, 
spherically symmetric and asymptotically flat spacetimes, to be 

CW(r) = A(r)/4 (74 = 71 6' (0/(74, (28) 

which coincides with the entropy at the horizon r = r^. 

In searching for generalizations of the C-function in higher curvature gravity, we continue 
to require that this function satisfies the usual properties [11], at least for the static and 
spherically symmetric specetimes : 

a) It can be evaluated on any spherical surface concentric with the horizon. 

b) When evaluated on the horizon of a black hole it reduces to the entropy. 

c) If certain physical (e.g. null or weak) and boundary (e.g. asymptotically flat or AdS) 
conditions are satisfied, then C is a non- decreasing function along the outward 
radial direction. 

Hence, let us take our proposed C-function of the pure Lovelock gravity of order m to be 
equal to the expression (24) evaluated on any spherical surface concentric with the horizon, 
i.e. we write 

^ 4nmn\ am ^n-2m+2 /X _ ^29) 

^ ' (n-2m + 2)! ^ ' ^ ' 

We have shown in the previous section that 6 is monotonically increasing function of r in 
pure Lovelock gravity of order m as long as am > and the matter content satisfies the 
null energy condition. Hence C^^^ satisfies the above three conditions and can serve as a 
candidate for a well defined C-function. We call these functions C-functions of the first kind 
However, in the next section we will show that this is not the only C-function one can 
define and yet satisfy the conditions stated above. 

3.2 C-functions of the second kind 

In this section we show that another class of well defined C-functions exists. Our proposed 
form is motivated by expression (27) after dropping out the subscript H allowing the cal- 
culations of this quantity on any sphere concentric with the horizon. Hence, our second 

''The C-function of the first kind was noted previously without proof in [28]. 
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C-function takes the form 



^ - (n-2m + 2)! 1, ¥{r) ) ' ^ ^ 

Taking constant time slices of the metric (5), we notice that the term a^6'^ is proportional 
to the extrinsic curvature of constant r surfaces. Taking the normal to be n = dr /a{r) we 
obtain ^ 

C("*) ~ 6" (Rf~ , (31) 

and 

C^"^) {R-K'' + KijK'^y^ , (32) 

where R and = Vifij are the intrinsic and extrinsic curvature of spheres. 

Obviously, the function (30) also gives the correct form of the entropy when evaluated 
on the horizon, — 0. Now we turn to the question whether C'^™') satisfies condition (c) 
above. In the following we show that this function, indeed, increases with radius provided 
that the matter content satisfies the weak energy condition. 

To make the notation compact, let us drop out the preceding numerical coefficients in 
(30) and consider instead the function 

F(r) = K (l - h'^y' , (33) 

where k — am'n\/2{n — 2m + 1)!. If this function were non-monotonic, then we could find a 
radius Tq with th < Tq < oo such that F'{ro) — 0. We start by writing the function F in the 
form F = b{b"rjn,n) / (1 ~ ct^^'^); where r^^„ is given by eq. (21). By direct calculations 
we find that 

dF ^ b'b^r^,n + b{b^rm,n)' 2b^+'r^,r,{aa'b'' + aH'b") 

dr ' ^ ' 

Further, we use eq. (20) and the t — t component of the equations of motion (10) to substitute 
for the quantities (6" r^^„)' and (a a' b'"^ + €? b' b"), respectively which yields 

dF b'b^ fn-m + 1 , m-1 ^ \ 

-rm,nH bp] . (35) 



dr 1 — a'^b'^ \ m ' m 

However, we have shown before that 6' > ( b is monotonic), and the weak energy condition 
was enough to prove Vm,n > 0. We conclude immediately that the r.h.s of (35) is positive 
definite, and hence there is no solution to dF/dr = 0. This proves that the functions C^"^\r) 
are C-functions for pure Lovelock gravity coupled to matter that satisfies the weak energy 
condition. 



3.3 Example: the C-functions in the vacuum solution of pure 
Lovelock gravity 

We can check our results by looking at Vacuum solutions of pure Lovelock gravity. These 
solutions can be obtained from eq. (22) by putting p — , i.e. the vacuum solution is given 
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by [24, 27] 

a\r) = l-(—) ,b{r)^r. (36) 



r 

The C-functions of the first and second kind for these spacetimes are given by 

Mm) ^ 47rmn!Q:A n-2m+2 

(n-2m + 2)! 

~, , , , \ (n-2m+2)(m-l)/m 

^ 4 TT m n! Q;^Q„ (n-2m+2)/m (n-2m+2)(m-l)/m /gyx 

(n-2m + 2)! " ' ^ ^ 

Prom the above equations we see that both C("^) and C^""^ are monotonic functions of the 
radial coordinate, and both reduce to the entropy when evaluated on the horizon. We also 
see that C^"*) < C^'") for all r > m- 



4 General Lovelock gravity, entropy and C-functions 

Our results in section 3 were for pure Lovelock theories, with the coefficient of only one of 
the Lovelock terms in the Lagrangian nonzero. Now we want to ask whether these results 
hold in a general Lovelock gravity theory. 

The Lagrangian for general Lovelock gravity is given by £ = J2m=f ctm^m, where Cm are 
given by eq. (4) , and we drop the cosmological constant as we are interested in asymptotically 
fiat solutions. The equations of motion read 

[(D-l)/2] 

GJ^ E OimGj^'^'^/2^+^^-Tj . (38) 

m=l 

Similarly the entropy is given by 5" = Z]|i=i^'*^^^ g(rn) ^ with S^"^^ given by eq. (24). 

In the previous section we proved that 6 is a monotonic function of r in all pure Lovelock 
theories. We can show that this result still holds for general Lovelock theories. As before, 
we take the combination 

[(D-l)/2] 

E (^t* - G/ ^'"^) /2™+^ = T; - Tl = TabC , (39) 

m=l 

and using eq. (10) we obtain 

a' b"/b'' = '^^^^ . (40) 

Z^m=l (n-2m+l)! \ b"^ J 

Assuming the positivity of the coupling constants am for all orders, and using the same 
reasoning as in pure Lovelock gravities, we conclude that b is also monotonic in general 
Lovelock gravity coupled to matter that satisfies the null energy condition. 
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4.1 C-functions in general Lovelock gravity 

For general Lovelock gravity, one can define either the general C-function of the first or 
second kind. In the first case the C-function is given by 



[(■D-l)/2] [(D-l)/2] 7,-2m/„A 

ii i'l (n-2m + 2)! 



Taking the derivative w.r.t r we obtain 



■^^ {n - 2m + l)\ 



(41) 



(42) 



As we showed before 6' 7^ for th < r < oo. Moreover, the positivity of the couphng 

constants am ensures that there is no solution to the polynomial under the sum. This proves 
the monotonicity of the general C-functions of the first kind provided that the matter content 
satisfies the null energy condition. 

In the same way we define the general C-function of the second kind to be 



[{D-l)/2] 
m=l 



(43) 



However, testing the monotonicity of this function is generally complicated. In the following 
we restrict our analysis to the case of Gauss- Bonnet gravity. 



4.2 C-function of the second kind in Gauss-Bonnet gravity 

The C-function of the second kind in D — n + 2 dimensional Gauss-Bonnet gravity reads 



Cgb =47rQ„ [Q;i6" + 2n(n- 1)q;2 6""^ (l - aH'^) 



(44) 



To prove the monotonicity of this function we proceed as we did before. We define the 
function 

F{r) = 6" + 2 n (n - 1) q; 6"-^ (l - b'^) , (45) 
where a — a2/ai. We then ask whether it is possible to find solutions to dF/dr — where 



^^nb'b^-^ 
dr 



b^ + 2{n- l)a ((n - 2) (l - b'^) - 2b (a^ b" + aa' b')) 



(46) 



Using the equations of motion (10) and (38), the relation a^Trr + Tu/a^ — T"-''^aCb: with 
being a null vector, and the definition of F above we find 



aH" + aa'b' 



[6(F-6")/4+ (n-2)(n-3)(F-6")V8n(n- 1)6"-^ -nlafe^+^Ttt/a^] 

n q; 6" (n - 2) a{F - b"") 



(47) 

Substituting this expression into eq.(46), we find the the solutions of F' = are given by 
the solutions of the equation 



6 6" 



(n- l)(n-2) 



F + 



(n + l)(n + 2) SnnlaTtt 
(n-l)(n-2) ^ {n-2)~^ 



6'" = , 



(48) 



10 



from which we see immediately that there are no real solutions for n > 3 if a > 0, and 
Ttt = P ^ 0, i.e. for matter content that satisfies the weak energy condition. This proves the 
monotonicity of the C-function of the second kind in Gauss-Bonnet gravity. 

Going to higher order Lovelock gravity, at least according to the present method, requires 
thorough analysis of higher degree polynomials. We will not attempt to carry this out here. 
Instead, in the next section we use numerical techniques to study the monotonicity of the 
C-function of the second kind. For a particular type of matter satisfying the weak energy 
condition, we will verify that C is monotonic for Gauss-Bonnet and see that it is also mono- 
tonic including the third order Lovelock term with positive coefficient. This result suggests 
that it might be possible to improve on the results in this section and show monotonicity of 
C{r) for all Lovelock gravity theories with coefficients > 0. 



4.3 Numerical example: Gauss-Bonnet and third order Lovelock 

gravity 

In this section we work out various numerical examples that show the monotonicity of C- 
functions of the second kind. In the following we consider general Lovelock gravity in D 
dimensions coupled to two abehan gauge fields with a — 1,2 and a massless scalar 
modulus field 0. This has been recentely studied in the context of the attractor mechanism 
[29, 30]. 

The modulus scalar has vanishing potential, but couples to the gauge field kinetic terms 
through the matrix function fab{<P)- The action is given by 



[D/2] 

'^m -^m ^ >^juT"^ Jab\^J-^^u 

m 



(49) 



where /i,i/ = 0, — 1. Using the static and spherically symmetric ansatz (5), we look for 
solutions to the resulting equations of motion (38), where their explicit form is given for h-D 
in a previous work [30]. The equations of motion for the gauge fields 9^ {y/—g fabF"'^") — 
may be solved by taking the field strengths to be of the form 

fab Q 

pa ^ J^dt A dr , (50) 

where are the electric charges, and the field dependent tensor f"'^{4)) is the inverse of the 
tensor coupling /a&(0) that appears in the Lagrangian. With this form for the field strength, 
the energy-momentum tensor for the gauge fields can be written in terms of the effective 
potential 

v,^ = r\(t>)QcQd, (51) 

and hence 

T/ = -a^(W-^,„, 
T; = a\dAf-^^. (52) 



2 /a M\2 ^eff 
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Also, the effective potential acts as a potential in the equation of motion for the modulus 
scalar, which is given by 



2 6« 



(53) 



This equation may be solved by a constant value of the scalar field if this value represents 
a critical point of the effective potential, i.e. Vgj^(0) = 0. Given that the scalar field is 
constant throughout the spacetime one obtains a simple solution to the equations of motion, 
namely summing over all Lovelock orders m in eq. (20) one obtains non-extremal black hole 
solution whose metric functions are solutions of the equation 



[{D-l)/2] , 2n-2m 

^ 2(n-2m + l)! 



m=l 



l+n—2m 



1 



n 



1- - 



r 



n— 1 



(54) 



and b{r) = r, where th is the outer horizon radius of the black hole. ^ 

Now we want to consider solutions with non-constant 0. To this end, we take small 
perturbations of the scalar field near the horizon where the metric functions are approximated 

by 

a^(r) p(r - th) , b{r)^rH, (55) 



and p is given by 



where 



P 



dr 



Si - Kff(0)/r5 



(56) 



Si = 



So = 



i^, 2(n-2m)!' 

2(n-2m+l)!- 



(57) 



Considering the scalar field perturbation 0(r) = -|- e0i(r), where e << 1, we find that 
the first order perturbative equation in the near horizon region is then given by 



(r - rH)^;' + 0; - 



2r^V 



0, 



(58) 



where = V^"(0). The well behaved solution for linearized perturbations of the scalar field 
is then given by 



/3 2{r-rn) 



(59) 



V P 

where Iq is the modified Bessel's function of the first kind and E is an integration constant. 

To this end, one can use the solution to the scalar field perturbation as initial condition to 
the full non-linear system. In the following we numerically integrate the non-linear equations 



^For explicit expression in Gauss-Bonnet gravity see [30]. 
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(38) and (53) using the Rung-Kutta method. We also take the couphngs of the scalar field 
to the gauge fields to be 

fabi^)=e-^^'^Sab, (60) 
from which we find immediately that the effective potential is given by 

Vesi<P) = e"'^Ql + e^'''Ql. (61) 




Figure 2: Numerical results for the C-functions of the first and second kind in D = 5 Gauss- 
Bonnet gravity. We choose rn = 1.26, coupling constants ai = a2 = 1, charges Qei = l/V^ 
and Qe2 = V^, 7i = —72 = 2.0 and 6r = 0.01 in our numerical scheme. We also take E = 0.4 
in the first two figures, while in the third figure we use E = 0.0 , 0.2 , 0.4 for the solid, dashed 
and dotted lines respectively. 

In figure (2) we compare the C functions of the first and second kind in Gauss-Bonnet 
gravity in five dimensions. We choose Th = 1.26, coupling constants ai = 02 = 1, charges 



13 
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Figure 3: Numerical results for the C-functions of the first and second kind in D = 7 third 
order Lovelock gravity. We choose th = 1.26, coupling constants ai = 0:2 = as = 1, charges 
Qei = I/V2 and Qe2 = v^, 71 = -72 = 2.0. We also take 6r = 0.01 and E = 0.4 in our 
numerical scheme. 

Qei = I/V2 and Qe2 = V2, 7i = -72 = 2.0. We denote the proximity to the horizon 
by the parameter 6r = (r^ — rH)/rj, this is where the initial conditions are set using the 
pertarbative results above, and we take it to be 0.01 in our numerical scheme. We also take 
E = 0.4 in the first two figures. Figure (a) shows the monotonic behavior of the C-functions 
of the first kind. On the contrary. Figure (b) shows the non-monotonic behavior of the term 
12 6(1 — 6'^) that appears in the C-function of the second kind in eq. (44). However, as it is 
clear form figure (c), when we add up the term the overall function restores its monotonic 
behavior. In addition. In figure (c) we compare the C-function for different values of the 
constant E, we take E = 0.0 , 0.2 , 0.4 for the solid, dashed and dotted fines respectively. We 
notice that the ultraviolet value of C decreases as we increase the value of the constant E, 
or in other words as we increase the asymptotic value of the scalar modulus. 

In figure (3) we show the results for the case of third order Lovelock gravity in seven 
dimensional spacetime. We see a similar behavior to the case oi D = 5: although the third 
term in C is decreasing, the overall function is monotonically increasing in r. 
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Higher order Lovelock terms have also been considered numerically up to the fifth order, 
and all results show monotonic behavior for C{r). 

5 Conclusion 

In this paper wc have constructed two different C-functions for the static, spherically sym- 
metric black holes in Lovelock gravity. This construction was inspired by Wald's expression 
for the entropy of stationary black holes applied to Lovelock gravity. Although this expres- 
sion is given as an integral over the induced metric on the Killing horizon, we were able to 
show that extending this expression in two different ways, by evaluating it on any spherical 
surface concentric with the horizon, gives the desired C-functions. These functions have 
different asymptotic values, but they degenerate at the horizon to the entropy of the black 
hole. 

In the case of pure Lovelock gravity of order m, the expression of the C-function of 
the first kind is simply proportional to 5""^™"+^, while that of the second kind is given by 
the former expression multiplied by the factor (1 — a^6'^)™^^ which was shown to contain 
contributions from the extrinsic curvatures of an n-sphere embedded in n + 1 dimensional 
space. We have also proven the monotonicity of these functions provided that the sapce is 
asymptotically fiat, and the matter content satisfies the null and the weak energy conditions 
for the first and second C-functions, respectively. It is worth mentioning that one can still 
show the monotonicity of the C-function of the first kind if we replace the asymptotically 
fiat by asymptotically AdS space since the later satisfies the null energy condition. 

In a general Lovelock gravity. It is natural to expect that the C-functions can be obtained 
by summing over pure C-functions of different orders. Although we proved the monotonicity 
of the general C-function of the first kind, we proved the monotonicity of the C-function of 
the second kind only in the case of Gauss-Bonnet gravity. It is still not obvious how to give a 
similar proof in the case of higher order gravity. However, we have given a numerical example 
in third order Lovelock gravity that indicates that the monotonicity of C{r) may still hold 
in general. We have also checked numerically that the results hold in similar examples in 
higher order Lovelock theories. It is worth noting that a quasi-local mass function in Gauss- 
Bonnet gravity was defined in [31] which also exhibits the monotonicity behavior under the 
dominant energy condition. 

The existence of two different C-functions raises the question whether there is some 
reason to prefer one over the other on physical grounds. The answer to this question relies 
on the existence of a covariant formulation of the C-function, which may reduce to one of 
the C-functions defined in this work, when evaluated in static and spherically symmetric 
spacetimes. In turn, the existence of such a covariant function would establish a second law 
of black hole mechanics (which was estabhshed for Einstein gravity in [32] ) in Lovelock 
gravity: if certain energy condition is satisfied, then the entropy of a black hole can never 
decrease. In other words, it may be that the second law of black hole mechanics selects the 
C-function that respects the law. Work along these lines is in progress. 
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